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In this paper we study updating effects on cellular automata rule space. We consider a subset of 
6144 order-3 automata from the space of 262144 bi-dimensional outer-totalistic rules. We compare 
synchronous to asynchronous and sequential updatings. Focusing on two automata, we discuss how 
update changes destroy typical structures of these rules. Besides, we show that the first-order phase- 
transition in the multiverse of synchronous cellular automata, revealed with the use of a recently 
introduced control parameter, seems to be robust not only to changes in update schema but also to 
different initial densities. 
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INTRODUCTION 

The cellular automata (CA) was created in the 1940s 
by Stanislaw Ulam and John von Neumann to study self- 
reproducing biological systems. Since then, deterministic 
CA gave origin to a precific field of study itself. In par¬ 
ticular, Wolfram [lj, Q has shown the whole variety of 
behaviors and patterns exhibited by automata rules. He 
argued that any CA rule could belong to one of the fol¬ 
lowing behavioral classes: fixed point (Class I), periodic 
(Class II), chaotic (Class III) and complex (Class IV). 
Besides, Langton [3] claimed the existence of a control 
parameter (A) that would classify CA rules in general 
and reveal a phase transition in CA rule space if plotted 
against a suitable but yet “unknown” order parameter. 

A recent attempt to determine the type of phase tran¬ 
sition CA rules are related to was made by Reia and 
Kinouchi [fj . We have proposed a control parameter a 
which classifies the 2 dimensional (2 D) automata rules 
in two regimes: a dead absorbing state (or vacuum) and 
a regime of high density of alive sites. The transition 
between both regimes was found to be of first-order with 
a metastable region, and the most interesting rules. We 
like to think, along other authors, that the rule space, 
with its diverse dynamical laws, represent a toy model 
for the multiverse jl|. In this sense, we found that com¬ 
plex “universes” like Game of Life [6| are special because 
they lie in the metastable region of the transition, which 
is a better description than the edge of chaos concept for 
CA 0. 

Cellular automata are not only used to study criti¬ 
cal properties of deterministic abstract dynamical sys¬ 
tems, but also to study natural and social phenomena. 
In general, CA can be applied where there are inter¬ 
est on dynamical evolution, global pattern formation 
and statistical properties of multiple interacting agents. 
This agent-base modelling spread over a wide range 
of science domains and is very important to computa¬ 
tional simulations of complex systems. In these cases, 
however, the automaton is usually probabilistic. We 
find, for instance, applications from ferromagnetism, like 


Ising 0| and Blume-Capel@ models, to social interac¬ 
tion modelling, like majority vote model 0 and Ax¬ 


elrod model 11], Examples of further applications of 
probabilistic CA can still be found in lj, |13| ■ 

The equivalence between computer simulations and 
real problems was addressed by Huberman and 
Glance 0 - These authors emphasize that one must be 
concerned about the effects of the update procedure on 
the global results of the simulation. According the au¬ 
thors, this concern is necessary to obtain a reliable in¬ 
sight into the real problem. Specifically, Huberman and 
Glance studied the Prisioner’s Dilemma in a square lat¬ 
tice with synchronous and asynchronous updating and 
showed the results strongly differ. 

In this perspective, this paper presents a study of the 
effects of updating schemas and different initial densities 
on the CA phase transition reported in 0. In particu¬ 
lar, we study how asynchronous and sequential updating 
affect the 2D CA rule space. Regarding the results re¬ 
ported in 0, in which authors investigate how Game 
of Life steady states are related to its initial density and 
asynchronity, we also study the stability of 2D patterns 
which have similar densities to the fixed point densities 
predicted by mean field analysis. Our results suggest that 
although nonsynchronous updating destroys characteris¬ 
tic structures found in synchronous updating, which is 


in agreement with [14|, [l6j,and that although the initial 


density can act in favor of dead or alive domains, the 
first-order phase-transition found in 0 seems to remain 
unchanged. 


RULES AND UPDATING SCHEMAS 

We will work in the rule space of binary outer-totalistic 
2D CA. The most known example of such rules is the 
Game of Life, where each site is either dead (0) or alive 
(1). From any initial configuration, the state of a given 
site is updated according its actual state s and the num¬ 
ber h of alive neighbors in its Moore neighborhood. Pre¬ 
cisely, a dead site becomes alive (“born B ”) only if it 
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has three alive neighbors, and an alive cell remains alive 
(“survive S ”) only if it has two or three alive neighbors - 
all other situations lead to a dead site at next time step. 
An usual code for this rule is B3S23. 

Our results are obtained in a 2D lattice with periodic 
boundary conditions. All results reported here are from 
the steady steady of the system, at t = 5000. The in¬ 
fluence of update schemas can be viewed in Fig. [Q In 
Fig. 1(a) we have the usual synchronous update which 
leads to the well known structures of Game of Life E3- 

The implication of an asynchronous update for Life is 
seen in Fig. |l(b)| The asynchronous update is realized 
by choosing a site randomly and updating it according 
to the Life’s rule. In this case, a time step is composed 
of N = L m x L n interactions and we see that the pattern 
is a labyrinth phase: the Life’s characteristic configura¬ 
tions found in synchronous update are destroyed and a 
different pattern is established. We notice the number 
of alive sites is higher than seen in Fig. 1(a) This sug¬ 
gests that asynchronous updating acts in favor of alive 
sites by allowing the formation of a stationary labyrinth 
phase with p^ = 0.436. A more detailed study about 
properties of Life under asynchronous updating can be 
found in [30. 

Fig. 1 1(c)] shows a snapshot of a pattern created by se¬ 
quential update, where sites are updated in an ordered 
way: sites of the first line of the lattice are firstly up¬ 
dated from left to right, then sites of the second line, and 
so forth. This updating schema seems to conduce the 
system to a similar labyrinth phase found in the asyn¬ 
chronous case but with a lower density, poo = 0.408. 
This lower number of alive sites reflects the existence 
of domains of dead sites that prevents the establishment 
of denser structures of alive sites. Curiously, the value 


Poo = 0.408 is very close to the site percolation threshold 
p c = 0.407 for a square lattice with Moore neighbor¬ 
ed m. 


MEAN-FIELD ANALYSIS 

Although we are studying bi-dimensional rules, it is 
interesting to consider if mean field analysis can give any 
insight for the 2D case. Each rule define an unique au¬ 
tomaton and is given by i?[s, h\. The rule R[s, h\ specifies 
the next state of a cell in the state s (s = 0 for a dead 
site, s = 1 for an alive site) with h alive neighbors. 

In this approach, we want to calculate the density of 
alive cells p(t) [17,|l9|. By neglecting spatial correlations, 
the time evolution expression for p(t) is given by 

1 8 

p(i + 1) = EE R[s,h]P t (s,h), (1) 

s=0 h —0 

where the term P t (s, h ) is the probability of a site s hav¬ 
ing h alive neighbors at time t. 

By using P t (s,h) = P(s)P(h,t ), it can be shown that 
Eq. (]T|) results in 


m (p) = Pit + 1) = (1 -Pi*)) ELo R[0,h]P(h,t) 

+pit ) (2) 

Here, P(h,t) is the probability of finding a number h of 
alive neighbors, which is assumed to follow a binomial 
distribution: P(/i,t) = C®p(£) ft, (l — p{t)) 8 ~ h . 

The Life’s rule is given by P[0,3] = 1, P[l,2] = 1, 
P[l, 3] = 1 and R[s, h] = 0 for other values of s, h. The 
application of this rule in Eq. (0 produces the map: 



FIG. 1: Snapshot at steady state of Life with syn¬ 
chronous 1(a) | asynchronous 1(b) and sequential 1(c) updates 
in a 2D lattice with L m x L n = 100 x 50. The system 
evolved to steady state from a randomly initial condition with 
p(0) = 0.5. 


M(p) = 28p(f) 3 (l - p(t)) 5 ( 3 - p{t)). (3) 

This map is a polynomial dominated by a order-3 term 
pit) 3 when p —> 0 (which was confirmed numerically in 
[19(| 4, and we call it a order-3 rule. This limit behavior 
is important to determine the stable points of the return 
map in Eq. ©. 

The map of Eq. © can be calculated for any rule and 
gives the temporal evolution of the density p{t) given 
pi 0). It is valid and describes very well CA rules embed¬ 
ded in random networks with eight neighbours. However, 
it fails in 2D for several important rules, like Life. This 
induced researchers to think that it is useless to under¬ 
stand complex rules. We have shown jj] that this is not 
the case: the mean-field fixed points can aid to define a 
control parameter a that locate complex rules inside the 
metastable region of a phase transition. Complex rules 
are complex because in 2D there is nucleation, compe¬ 
tition and coexistence of domains (formed by the dense 
fixed point phase) with the zero phase. 
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RESULTS 

In this context, we focus our attention only on the 
order-3 rules from Eq. 0 because they are the most 
dynamically interesting. These rules generally have two 
stable fixed points: p = 0 and p = p*. Reia and Ki- 
nouchi 3 have claimed that 2D system presents alive 
and dead domains with similar densities as the mean- 
field fixed points. Because of that, we study the subset 
of 6144 (including Life) order-3 automata in the space 
of 262144 outer-totalistic CA rules to analyze the coex¬ 
istence and competition between these domains. From 
this subset, 2643 rules were excluded from our analysis 
by not having a stable fixed point p*. Thus, our sample 
is composed of 3501 rules. 

In order to see if the square lattice regards any resem¬ 
blance with mean-field predictions, we plot P 2 d/P* at the 
stationary state as function of a control parameter a in 
Fig. [2j where p 2 D is density of alive sites in the lattice. 
The parameter a was devised to describe the growth rate 
of the vacuum phase at the interface of domains. Con¬ 
sidering a square lattice composed of alive bubbles with 
density p = p* and dead bubbles with p = 0, one can 
show that density at interface between both bubbles is 
p*/ 2. Hence, the expansion and growing of dead sites 
can be inferred by defining the zero phase growth rate at 
interface Q: 


1 - M(p*/2) 

l-p*/2 


(4) 


The zero phase expands for er > 1 and contracts for 
a < 1. The critical value is a c = 1 and is at the tran¬ 
sition region of Fig. [2] Fig. [2] shows a large number of 
rules where 2D correlations are weak, which results in 
P 2 D = p*■ For the other rules, we distinguish three cases: 
(I) spatial correlations promote overactivy, resulting in 
P 2 d/p* > 1; (II) the mean-field prediction overestimates 
the density of alive cells, which results in a coexistence 
and competition between alive and dead domains with 
0 < P 2 d/pmf < 1; (HI) despite the stable p* predicted 
by mean-field calculations, bubbles of zero phase expand 
over the entire lattice, resulting in an absorbing phase 
with p 2 D = 0. 

What we also see in Fig. [2] is that different updating 
schemas do not affect the phase transition as it does to in¬ 
dividual 2D structural patterns. In fact, we observe that 
the destruction of structural patterns reduces the number 
of rules exhibiting coexistence and competition between 
domains (case(II)). As a consequence, although the tran¬ 
sition region is less populated in asynchronous |2(b)| and 
sequential 2(c) than synchronous 2(a) update, the dead 
absorbing state and the high density state regimes can 
be identified by using a. 

We have mentioned that the order-3 rules generally 
have two stable fixed points. The results shown in Fig. [2] 





FIG. 2: Relation between p 2 D/p* and a. Results were ob¬ 
tained for L = 100 and averaged over 30 runs. We chose 
p( 0) = p* to see the effect of a on competition of alive and 
dead cells in 2D with synchronous update in panel 2(a) with 
asynchronous update in panel 2(b)[ and sequential update in 
panel 2(c) The Life a = 1.006 and its position is marked with 
a large grey dot in each panel (pointed by the grey arrow). 
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report the phase transition found for p(0) = p*. In Fig.[3j 
we show the dependence of p 2 D / P* on cr when the initial 
density is p(0) = 0 + A p, where A p = 0.05, under the 
three update schemas. We observe that the transition 
region is broadened and shifted to left in all schemas. We 
can note that the three cases already discussed are also 
seen here. However, we also observe that the low initial 
density reduces the number of rules having P 2 d/p* > 1 
(case (I)), mainly in synchronous (3(a)) and sequential 
(3(c)) updates. This occurs because the high number 
of bubbles of zero phase found at the beginning of the 
simulation remains stable during the dynamics. Then, 
the growth of alive sites is constrained and there is a 
lower number or rules exhibiting overactivity. 

Considering that the low initial density acts in favor 
of bubbles of zero phase, one could expect that the dead 
sites would prevail over alive domains in a scenario of co¬ 
existence and competition. This expectation is confirmed 
by the higher number of rules having an absorbing phase 
P 2 D = 0 (case(III)) in Fig. [3] In fact, it also explains 
the higher number of rules in region 0 < P 2 d/P* < 1 
(case(II)). Specifically, the initial density is too low that 
even when a < 1 the rules are attracted to the stable 
fixed point p = 0. As a consequence, dead domains be¬ 
come stable and the value of P 2 d/ P* of these rules de¬ 
crease. 

As a matter of illustration, we show in Fig.0]2£) steady 
states for a rule with a < 1 whose value P 2 d/P* ~ 1 decay 
when the initial density is decreased. This rule is defined 
by H4678S'2345, which has cr = 0.956 and two stable 
fixed points: p — 0 and p* = 0.615. We see in Figs. |4(a)| 
|4(b)| and 4(c) the steady state with synchronous, asyn¬ 
chronous and sequential updates respectively. Note that 
the initial condition p(0) = p* allows a strong coexis¬ 
tence and competition between alive and dead domains 
in a synchronous updating, resulting in p^ = 0.491. As 
already discussed for Life, the asynchronous and sequen¬ 
tial updates acts in favor of alive sites. Consequently, 
alive sites prevail over dead domains and results in a sta¬ 
tionary density p^ = 0.678 in both cases. 

We see in Figs. 4(d) 4(e) and |4(f)| the steady states for 
the three updates when p(0) = 0 + Ap, where A p = 0.05. 
We can note that the high number of dead domains found 
at the beginning of the dynamics expands over the entire 
lattice. This results in a small number of nuclei of alive 
sites, which results in p^ ss 0 in the three schemas. This 
rule illustrates the underlying properties of the a < 1 
rules that broaden and shift the transition region to left 
when the stability of the absorbing steady state is studied 
in Fig. [3] 

The steady states of the system when an intermediary 
value p(0) = 0.200 is assigned to the initial density is 
seen in Figs. 4(g) 4(h) and 4(i)| under the three update 
schemas respectively. These snapshots depicts not only 
the strong coexistence and competition between dead and 
alive domains, which can be disturbed just by increasing 



(a) 




a 


(c) 


FIG. 3: Relation between P 2 d/p* and cr. Results were ob¬ 
tained for L = 100 and averaged over 30 runs. We chose 
p( 0) = 0 + A p, where A p = 0.05, to see the effect of a on 
competition of alive and dead cells in 2D with synchronous 
update in panel 3(a) with asynchronous update in panel|3(b)| 
and sequential update in panel [3(c) The Life cr = 1.006 and 
its position is marked with a large grey dot in each panel 
(pointed by the grey arrow). 
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or decreasing the initial density, but also the effect of the 
update in this subtle equilibrium. In Fig. |4(g)~| we observe 
the formation of stable structures like “crystals” (or “di¬ 
amonds”) of alive sites that persist until the end of the 
synchronous dynamics, with p^ = 0.170. In Fig. 4(h) 
the asynchronous update destroyed these typical struc¬ 
tures and favored the spread of alive domains, resulting 
in poo = 0.683. Curiously, Fig |4(i)| shows that sequential 
update promotes the establishment of stable structures 
of alive sites. However, in contrast to synchronous case, 
the lattice is poorly populated by alive domains, allow¬ 
ing the presence of greater regions of dead sites, but these 
nuclei of alive sites are denser than in synchronous case, 
which results in p^ = 0.140. 

It is noteworthy the results from Fates [THJ. In this 
work, the author studies the dependence of the steady 
state density poo on the initial density p(0) and on the 
degree of desynchronisation of the update. He finds that 
Poo of Life is not related to p(0) when the dynamical up¬ 
date is asynchronous. Our results suggest this is not the 
case for all automata rules. The snapshots at the steady 
state of the rule B4678S2345 showed in Fig. |4(b)| and 
Fig. 4(e) reveals that the initial density affects the final 
pattern in the lattice and the final density. This may hap¬ 
pen to the asynchronous rules that present a decay from 
P 2 d/p* ~ 1 in phase transition (Fig. |2(b)[ ) when the ini¬ 
tial density is changed (Fig. 3(b)), as already discussed. 
Despite that, the general aspects of phase transition re¬ 
mains qualitatively the same. 


CONCLUSION 


We study the effects of asynchronous and sequential 
updates in the dynamics of 3501 CA in a square lattice. 
By using Life as example, we see that updates changes 
destroys characteristic structures of this automaton and 
are favorable to the spread of alive sites. We have shown 
that the first-order phase-transition found with the use 
of the recently proposed control parameter a is observed 
not only in the synchronous but also in asynchronous and 
sequential updates. 

The parameter a measures the growth rate of the 
vacuum phase at the domains interface. Both updates 
schema also present a transition from a high density state 
to a dead absorbing state around cr c = 1, which indi¬ 
cates that the control parameter a, although heuristic, 
seems to be robust. Indeed, this transition was found 
to be qualitatively the same when the initial density was 
changed. The fact that the transition is robust and that 
complex rules lie near <t c = 1 illustrates the idea that, 
in the multiverse of CA, complex universes lie inside a 
phase transition region. 



FIG. 4: Snapshot at steady state of rule B4678S2345 with 
synchronous |4(a)| asynchronous 4(b) and sequential 4(c) up¬ 
dates in a 2D lattice with L m x L n = 100 x 50. The sys¬ 
tem evolved to steady state from a randomly initial condi¬ 
tion with p(0) = p*. Figs. 4(d)| 4(e) and 4(f) show steady 
states for the three schemas when p(0) = 0 + Ap, where 
A p = 0.05. Figs. 4(g) | 4(h) and 4(i) show steady states for up¬ 
dates synchronous, asynchronous and sequential, respectively, 
when p(0) = 0.200. 
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